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ABSTRACT 

Flares in Sagittarius A* are produced by hot plasmas within a few 
Schwarzschild radii of the supermassive black hole at the Galactic center. The 
recent detection of a correlation between the spectral index and flux during a 
near infrared (NIR) flare provides a means to conduct detailed investigations 
of the plasma heating and radiation processes. We study the evolution of the 
electron distribution function under the influence of a turbulent magnetic field 
in a hot collisionless plasma. The magnetic field, presumably generated through 
instabilities in the accretion flow, can both heat the plasma via resonant wave- 
particle coupling and cool the electrons via radiation. The electron distribution 
can generally be approximated as relativistic Maxwellian. To account for the 
observed correlation, we find that the magnetic field needs to be anti-correlated 
with the electron "temperature". NIR and X-ray light curves are produced for 
a cooling and a heating phase. The model predicts simultaneous flare activity 
in the NIR and X-ray bands, which can be compared with observations. These 
results can be applied to MHD simulations to study the radiative characteristics 
of collisionless plasmas, especially accretion flows in low-luminosity AGNs. 

Subject headings: acceleration of particles — black hole physics — Galaxy: center 
— plasmas — radiation mechanisms: thermal — turbulence 

1. Introduction 

It is generally accepted that the near-infrared (NIR) and X-ray flares in Sagittarius 
A*, the compact radio source associated with a mass M ~ 3 — 4 x 10 6 M supermassive 
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black hole at the Galactic Center (Schodel et al. 2002; Ghez et al. 2004; the Schwarzschild 
radius r s = 10 12 (M/3.4 x 10 6 M ) cm), are produced by relativistic electrons within a few 
Schwarzschild radii of the black hole via synchrotron and synchrotron self-Comptonization 
(SSC) processes, respectively (Baganoff et al. 2001, 2003; Genzel et al. 2003; Belanger 
et al. 2006; Eckart et al. 2004, 2006; Yusef-Zadeh et al. 2006a). Studying the details 
of the electron acceleration therefore plays an important role in revealing the properties 
of the flaring plasma and its interaction with the surroundings and the space-time created 
by the black hole, which may lead to a measurement of the spin with comprehensive flare 
observations and general relativistic MHD and ray-tracing simulations of the black hole 
accretion and the corresponding radiation transfer. 

The electrons are likely accelerated by turbulent plasma waves generated in an accretion 
torus via instabilities, which induce the release of gravitational energy of the plasma (Balbus 
& Hawley 1991; Tagger & Melia 2006). In the simplest case, in which there is no large scale 
magnetic field, the turbulent magnetic field B can play the dual role of both heating the 
electrons (presumably via resonant wave particle interactions) and cooling them (via SSC 
radiation), since the heating and cooling rates have different energy dependencies. Earlier, 
we showed that the steady-state electron spectrum can be approximated as a relativistic 
Maxwellian distribution with the "temperature" r y c m e c 2 depending on the plasma density n 
and the coherent length of the magnetic field, which should be comparable to the size of the 
flare region R. Simultaneous spectroscopic observations in the NIR and X-ray bands (with 
four measured quantities: fluxes and spectral indexes in both bands) can then be used to 
determine the four basic model parameters, namely B, R, r y c m e c 2 and n, and to test the 
model (Liu et al. 2006a, 2006b). 

These investigations, however, also indicated that the steady-state solutions are not al- 
ways applicable, especially for NIR flares with very soft spectra. A time-dependent approach 
to the evolution of the electron distribution under the influence of the turbulent magnetic 
field is required to study the flares in quantitative details. In light of the recent discovery of 
a correlation between the flux F v oc v a ~ l and the spectral index a in the NIR band (Ghez 
et al. 2005; Gillessen et al. 2006; Krabbe et al. 2006; Hornstein et al. 2006), we carry 
out such an investigation. We show that relativistic Maxwellian distributions in general give 
a good description of the electron spectra with the corresponding "temperature" evolution 
determined by the initial temperature, the magnetic field and its coherent length and the 
gas density. It is therefore possible to couple these results with MHD simulations to study 
the radiative characteristics of collisionless plasmas self-consistently. 

Several scenarios have been suggested to explain the observed correlation between flux 
and spectral index (Gillessen et al. 2006). We find here that the dynamical processes of 
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plasma heating and cooling are the most likely mechanism and the observation requires that 
the magnetic field be anti-correlated with the electron temperature. Since the evolution of the 
magnetic field depends on the turbulence generation mechanism, which needs to be addressed 
via MHD simulations, we calculate the evolution of the electron distribution function and 
the corresponding radiation spectrum for a cooling and a heating phase, where the magnetic 
field is set to be inversely proportional to the mean electron energy. The model naturally 
recovers the observed correlation between the NIR flux and spectral index and predicts the 
X-ray emission characteristics accompanying the NIR flares. Simultaneous flare observations 
in the NIR and X-ray bands can readily test the model and uncover the underlying physical 
processes producing these flares. 

The kinetic equation of the electron acceleration and its time-dependent solutions are 
described in § 2. In § 3 the model is applied to flares in Sagittarius A* by calculating the 
radiation from a flaring plasma with prescribed magnetic field evolutions. The implication 
and limitation of these results are discussed in § 4. 



2. Time Dependent Solutions of Electron Acceleration by Plasma Waves 

The theory of electron acceleration by a turbulent magnetic field in a hot collisionless 
plasma has been discussed by Liu et al. (2006b). Regardless the details of resonant wave- 
particle coupling, it was shown that the evolution of the electron distribution function N(j, t) 
can be described by the following kinetic equation 

ON _ d 
dt c?7 

where the synchrotron cooling and acceleration times are given, respectively, by 

r sy „( 7 ) = r / 7 = 9m^c 5 /4e 4 S 2 7 , (2) 
r ac = 2 7 2 / < A 7 A 7 /At >= d3Rc/v 2 A = 12irCinm p cR/ B 2 , (3) 

the Alfven velocity Va = -B/(47rnm p ) 1 / 2 , C\ is a dimensionless quantity of order 1, and we 
have assumed that there is no particle escape and therefore no return current in association 
with the electron acceleration 5 . The size of the flaring region, the magnetic field and gas 
density are given respectively by R, B and n. The scattering mean free path of the electrons, 



5 These processes are not expected to affect the results significantly as far as the escape time scale is much 
longer than r ac (Liu et al. 2006a). 
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C\R, is comparable to the source size. m e , m p , c, e, and 7 are the electron mass, proton mass, 
speed of light, elemental charge unit and the Lorentz factor of the electron, respectively. 

To demonstrate the behavior of the time-dependent solutions, we consider the simplest 
case, where r ac and r remain constant in time, and normalize the distribution function. Then 
the steady state is ^(7) = (^ 2 /2^) exp(— / y/ / y c ) , with 7 C = r /r ac = 3mgC 4 /167rCie 4 m p n R. 
Here the integral over 7 has been extended from 1 to to simplify the expression. We note 
that 7 C is independent of B because both the synchrotron cooling and acceleration rates 
are proportional to B 2 , while R and n control 7 C by affecting the turbulence heat rate via 
modification of the coherent length and Alfven velocity va, respectively. Once one specifies 
the initial electron distribution with r ac as the time unit, the solution only depends on r or 
the steady-state temperature j c . 

In what follows, the energy unit is assumed to be m e c 2 . Figure 1 shows the evolution 
of iV(7, t) starting with a relativistic Maxwellian distribution with an initial temperature 7* 
until t = 4r ac . The left panel has 7$ = 10 and 7 C = 200, which corresponds to a heating phase 
of the plasma. We note that the particle distributions are already close to the steady-state 
solution after the time step t = r ac . The right panel has 7$ = 200 and 7 C = 10 corresponding 
to a cooling phase. Because the synchrotron cooling (with a timescale 7 c T a c/7) dominates 
above 7 C , the particle distribution evolves faster than in the heating phase (in unit of r ac ). We 
also note the pileup of electrons to a nearly monotonic distribution in the early section of this 
cooling phase, due to more rapid cooling of higher energy electrons. Later the distribution 
is broadened by the diffusion term in the kinetic equation (the first term on the right hand 
side of eq. [1]). 

Figure 2 shows the correlations between iV(7, t) and the electron spectral index f3 = 
<91niV(7,t)/<91n7 at 7 = 50 (left) and 100 (right) for the two runs in Figure 1. Here the 
solid and dashed lines are for the heating and cooling phases, respectively (one may consider 
the time t as a parametric variable). The rising of (3 to above 2 is due to the pileup of 
electrons in the early cooling phase. The correlations are quite different for the two phases 
and energies. For a given magnetic field, these correlations mimic similar correlations in the 
synchrotron radiation, which we investigate in § 3. 

To characterize the evolution of N(j, t) in general, we next consider the evolution of 
the mean energy < 7 > (t) — J r yN( r y,t)d r y/ J iV(7, t)^. Because the heating of electrons 
by turbulence depends on the derivative of N(j,t) with respect to 7 via the diffusion term, 
the evolution of < 7 > can be very complicated. Fortunately for any smooth function Nfr) 
one may approximate this derivative term with a heating term, i.e. the total energy change 
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rate 

< 7 > 2 (2 + a) < 7 > 
< 7 > ~ ! + ^ i ! — , (4) 

7~ac7c ^ac 

where a " ' " indicates a derivative with respect to time and a is a parameter to be determined 
by the steady state solution. Then we have 

< >t= (2 + a) 7c 7^ 

* 7» - (7» - 7c) exp[-(2 + a)*/ r ac ] ' 

In the steady state, < 7 > t = 37 c , therefore a — 1. 

However the above approximation is not good enough for the intermediate steps of the 
heating phase when the particles have a broad distribution (Fig. 1) and the heating by the 
diffusion term is more efficient. To fit the numerical results we set the "a" in the numerator 
of equation (5) equal to one to be consistent with the steady state solution and leave the 
"a" in the denominator as a free parameter: 



7 * 7i - (7i - 7c) exp[-(2 + a)t/r ac ] ' 1 ' 

We find that a = 1.6 and 1.0 give good fits to the < 7 > evolutions of the heating and 
cooling phases, respectively. Figure 3 shows the evolution of < 7 > and the corresponding 
fits for the heating (left) and cooling (right) phases. Figure 4 gives the ratio of < 7 > to the 
fit values for several initial and final temperatures. The relative error is within 20% for the 
heating phase and within 10% for the cooling. The error also increases with the increase of 
the dynamical range, which is the ratio of the initial temperature to steady state temperature 
for the cooling phase and vice versa for the heating phase. Equation (6) therefore gives a 
good description of the evolution of < 7 > under the influence of a turbulent magnetic field. 

We note that equation (6) can be generalized to address the heating and cooling of 
collisionless plasma by turbulent plasma waves. The heating time r ac (see eq. [3]) usually 
is a function of time. Taking into account the effects discussed in the previous paragraph, 
equation (5) leads to 

< 7 >t = *W* (7) 

MHD simulations can give the time evolution of B, its coherent length, and n. Equation (7) 
can then be used in these simulations to address the heating and cooling of electrons. In § 
3 we study the SSC emission of these electrons during flares in Sagittarius A*. 
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3. Synchrotron Emission and SSC 

We consider the case where the particle distribution is isotropic with respect to the 
turbulent magnetic field. Then the synchrotron flux density and emission coefficient at 
frequency v are given, respectively, by (Pacholczyk 1970) 

FM = £ v (u) = -^-^Bnj dry J M±-» 2 Y /2 N{l)F{x), (8) 

where 

v __ Aixm e c v 
X = V c = 3eB(l-/i 2 )V2 7 2 ' 

poo 

F(x) = x K 5/3 (z)dz~2.U95x 1 / 3 e~ 2x + 1.348x 1 / 2 e- x -°- 249x ~ 163 , (10) 

J X 

D, fx, and K 5 / 3 are the distance to the Galactic Center, the cosine of the angle between the 
magnetic field and line of sight, and the corresponding Bessel function, respectively, and the 
approximation for F(x) is accurate within 8%. We therefore have the spectral index in a 
given narrow frequency range a = d\n(i/F u )/d\ni/. 

As shown by Liu et al. (2006b), most of the synchrotron radiation is emitted in the 
optically thin region for flares in Sagittarius A*. For a uniform spherically symmetric source, 
the self-Comptonization flux density is then given by (Blumenthal & Gould 1970) 

Similarly one can define the X-ray spectral index ax = dln(z/Fx)/dlni/ . 

3.1. Flux and Spectral Index Correlations with a Constant Magnetic Field 

To demonstrate how the heating and cooling processes affect the evolution of the radi- 
ation spectrum, we first consider the cases where the magnetic field B and the total number 
of electrons 47rni? 3 /3 remain constant. Figure 5 shows the evolution of the normalized syn- 
chrotron flux density spectrum 

^ _ F u (u,t)V3m e c 2 D 2 f° f 1 21/2 
e v (v, t) = e 3 BnR s = J d ^ J d ^ -») t > F ^ x > ( 12 ) 

for the two runs in Figure 1 with B = 100 G. Compared to the evolution of N("f), the 
evolution of e v {y) is less dramatic because the synchrotron spectrum is dominated by emission 
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from the more energetic electrons, and the radiation spectra are very similar to thermal 
synchrotron spectra. However, the difference between the heating and cooling phases is 
obvious. 

Figure 6 shows the correlations between a and e„ at v = 1.4 x 10 14 Hz (left) and 
1.4 x 10 13 Hz (right). The dotted lines indicate the observed correlations with different back- 
ground subtraction methods (Gillessen et al. 2006). The cooling phase fits the observations 
marginally, while the heating phase predicts a correlation much weaker than is observed. 
We note that for constant B and nR 3 the observed flux density is proportional to e u . The 
correlation between e u and a can be compared with observations directly by adjusting the 
normalization factor BnR 3 . We also note that for given iV(7, t) the correlation between e„ 
and a only depends on u/B. Therefore for the same runs the correlations at 1.4 x 10 14 Hz 
with B = 1000 G will be the same as those shown in the right panel. 

In § 2 we showed that the evolution of iV(7, t) is very similar in the heating and cooling 
phases for different initial and final temperatures, and that adjusting these temperatures 
changes the starting and ending points of the correlation in the e u -a plane. The shape of 
the correlation, however, does not change dramatically. We therefore conclude that for a 
constant (and uniform) magnetic field the observed correlation between flux and spectral 
index in the NIR band can marginally be attributed to the cooling of a plasma heated up 
instantaneously at the onset of the flare. 

Adiabatic expansion or compression and Doppler boosting have also been suggested to 
explain the correlations between spectral index and flux during flares (Yusef-Zadeh et al. 
2006b; Gillessen et al. 2006). To produce NIR flares with very soft spectra the electrons 
producing the NIR emission need to have a sharp high energy cutoff. As shown by Liu et 
al. (2006a, 2006b) and discussed in § 2, a relativistic Maxwellian distribution gives the most 
natural approximation to the electron spectrum. For a temperature of 70 (Liu et al. 2006b) 
we have, 

e u (v) = 0.5x M I(x M ) , (13) 



where 



I(x M ) = 4.0505x M 1/6 (l + 0.40a; M 1/4 + 0.5316a; M 1/2 )exp(-1.8899a;^ 3 ), (14) 
xm = u/uq = 4:7im e ch , /3eB'jQ , (15) 



and 



1/3 0.1000^ 4 + 0.2658 
a = 1.833 - 0.6300^ - — ^ ^— T . 

xYf + 0.4000x^ 4 + 0.5316 



(16) 



The left panel of Figure 7 shows the correlation between e v and a. As mentioned above, 
when the magnetic field is constant, this can be compared with observations directly. Such 
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a correlation can be produced when a plasma is expanding or compressed adiabatically in 
a uniform magnetic field so that the temperature 70 changes. The theoretical prediction is 
clearly inconsistent with the observed results. A very unusual electron spectrum is needed to 
make the model of adiabatic expansion in a uniform magnetic field in line with observations. 

Because £ v /u 2 is a Lorentz invariant, F u /u 2 is approximately a constant under the 
Lorentz transform. The correlation between F v and a due to Doppler boosting effects is 
therefore given by the correlation between e u /x 2 M and a, which is much flatter than the 
correlation between e u and a and therefore can not explain the observations. 

We conclude that in a constant and uniform magnetic field only SSC cooling of a hot 
plasma heated instantaneously is marginally in agreement with the observed spectral index 
and flux correlation in the NIR band. Turbulent heating, adiabatic processes and Doppler 
effects produce a correlation which is much flatter than the observed results. 

3.2. Flux and Spectral Index Correlations with a Variable Magnetic Field 

To improve the model fitting to the observed correlations, we next consider the scenario 
where the magnetic field is variable. The correlation between flux and spectral index is 
depicted by the correlation between Be u and a. For a relativistic Maxwellian distribution 
with a temperature 70, if -B70 is a constant in time, Be v is proportional to xm^v As shown 
in the right panel of Figure 7, the correlation of the solid line gives a better fit to the 
observations than all models discussed above. We note, however, that for a increasing from 
-4 to 1, the temperature of the plasma has to increase by about three orders of magnitude 
(Liu et al. 2006b). 

In the more general cases where B'Jq (with p as a constant) does not change with time, 
Be v is proportional to xFj} 2 ^ 9, ] e v . For p = 0, we recover the result with a constant magnetic 
field. It is clear that p has to be positive and less than 2 to make the flux and spectral index 
correlation steeper, implying an anti-correlation between the magnetic field and the mean 
electron energy during the flare evolution. The dashed line in Figure 7 (right) corresponds 
to p = 4/3, which fits the correlation in the soft- and dim- state but predicts a spectrum, 
which is softer than the observed ones when the flux is high. On the other hand, the recent 
re-analyses of the Keck observations of flares from Sagittarius A* appear to favor such a 
scenario (Krabbe et al. 2006). 

To take into account the heating and cooling effects, we choose a magnetic field inversely 
proportional to < 7 > during the evolution of the flare, which makes the heating and cooling 
rates depend on the electron distribution. Compared with the heating and cooling processes 
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studied in § 2, one more parameter is needed to obtain the solutions: the initial magnetic 
field. This parameter also sets the characteristic timescale of the system, which is related to 
the initial value of r ac . Using an explicit number scheme with the time steps scaled as r ac (t), 
we calculate the evolution of the synchrotron spectrum for a cooling and heating phase. 

The dashed line in Figure 8 (left) shows the correlation between e u B and a for the cooling 
phase at v — 1.4 x 10 14 Hz, which is consistent with the observations. The lower magnetic 
field, the initial and final temperatures, are 4 Gauss, 2000 and 20, respectively. Then we 
have nR = 2.1 x 10 19 C 1 _1 cm~ 2 and r ac = 6.5 x 10 3 (-B/10 G) -2 mins. Here a relatively large 
temperature change is chosen to produce significant variation of the synchrotron spectral 
index. The dashed line in the right panel gives the evolution of the magnetic field. Although 
the initial magnetic field is relatively low, the steady state magnetic field is 400 Gauss, which 
is much higher than the typical value for the quiescent state. 

The dashed lines in Figure 9 show the light curves of the flux (left) and spectral index 
(right) for the cooling phase. We note that the correlation between the flux and spectral 
index mostly occurs between 130 to 160 minutes. Before 130 minutes, the system evolves 
slowly due to the relatively lower magnetic field (< 20 Gauss). Because the magnetic field 
in the accretion torus in Sagittarius A* is about a few tens of Gauss in the quiescent state, 
the evolution before 100 minutes is likely irrelevant to flares in Sagittarius A* since the 
magnetic field is below 10 Gauss during this period. From 100 to 170 minutes the model 
predicted correlation is quite in line with observations. The model not only produces the 
correct correlation slope, but also recovers the typical variation timescale of a few tens of 
minutes. After 170 minutes, the system reaches a steady state. The dips near 160 mins in 
the flux and spectral index light curves are related to the turn around of the spectral index 
and flux correlation before reaching the steady-state. These are caused by the sharp cutoff 
(sharper than an exponential cutoff) of the electron spectra in the cooling phase, which 
makes the spectra softer than that of a relativistic Maxwellian distribution. To produce a 
flare state with a = —0.5 and F u = 5 mJy, the total number of electrons involved in the 
flare 4nnR 3 /3 = 2.6 x 10 41 , giving rise to R = 5.5 x lO 10 ^ 172 cm and n = 3.7 x 10 8 C7 3/2 
cm -3 . These values agree with results of previous studies (Liu et al. 2006a, 2006b). 

The heating phase fits the observed correlation marginally (the solid line in the left panel 
of Figure 8). The initial magnetic field, the initial and final temperatures are 200 Gauss, 10 
and 1000, respectively, giving rise to nR = 4.1 x 10 17 Cf 1 cm~ 2 and r ac = 1.3 x 10 2 (£?/10 G)~ 2 
mins. The solid lines in the right panel of Figure 8 and in Figure 9 depict the evolution of 
the magnetic field and the flux and spectral index, respectively. The heating phase proceeds 
much faster than the cooling phase due to the relatively high initial magnetic field and final 
temperature. The flux reaches the peak value within ~ 6 minutes, after which the spectral 



-10- 



index saturates near 0.2, and both the magnetic field and flux decrease slowly. To produce 
even harder spectra, one has to increase the value of £> 7 2 at the steady state dramatically (Liu 
et al. 2006b), making the heating phase proceed on an even shorter timescale (oc 1/-B 2 7c ). 
A very weak initial magnetic field (< 1 Gauss) and a very high temperature (>> 1000) are 
required to produce hard NIR spectra and a flux rising time of a few minutes. To produce 
a flare state with a = —0.5 and F v = 5 mJy, the total number of electrons involved in the 
flare 47rni? 3 /3 = 4.7 x 10 42 , which is more than ten times larger than that in the cooling 
phase. This leads to R = 1.6 x 10 12 Cj /2 cm and n = 2.5 x 10 5 C~ 3/2 cm" 3 . 

In general the spectral index a of the above models is always less than 4/3. To produce 
even harder spectra, one has to introduce self-absorption effects. For an electron temperature 
7o — 200, the size of the emission region is then given by 

R ~ (F u /2-K lQ m e ) ll2 D/v = 3.7xl0 7 ( J D/8kpc)(F ;/ /5 mJy) 1/2 ( 7o /200)- 1/2 (z//1.4xl0 14 Hz)- 1 cm , 

which is much smaller than any relevant length scales. We therefore do not expect NIR flares 
with spectral indexes larger than 4/3. Any solid detection of very hard NIR spectra will rule 
out this model and may be in conflict with any synchrotron models. 



3.3. Flux and Spectral Index Correlation in the X-ray Band 

The spectrum of the Comptonization component and its evolution should be quite dif- 
ferent for the two phases, and may be used to distinguish between them. Figure 10 shows 
the evolution of the radiation spectrum during the cooling (left) and heating (right) phases. 
Besides indicating the initial and steady-state spectra, we highlight the spectral evolution 
when the correlation between NIR flux and spectral index are prominent, corresponding 
to the shaded periods in Figure 9 and Figure 11, where the spectra are indicated by filled 
circles. For the model parameters chosen in §3.2, because the electron column depth nR in 
the cooling phase is 50 times that in the heating phase while their synchrotron luminosities 
are comparable, much more high energy emission is produced via SSC in the cooling phase. 

In analogy to equation (12), one can define a dimensionless SSC flux density spectrum 

, , 6V3myD 2 F x (u) f°° Nh) f°° , ,eJi/) ( n , v t 2 . v 2 \ 

ne 7 n 2 R 4 B J T J " v V 47V 27V / 

(17) 

The left panel of Figure 11 shows the evolution of Be x at 0.5 x 10 18 (thin lines) and 2 x 
10 18 Hz (thick lines) for the cooling (dashed lines) and heating (solid lines) phases. As 
expected, higher energy emission precedes (lags) lower energy emission slightly during the 
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cooling (heating) phase because they are produced by more energetic electrons. However 
the difference in the peak times is within 10 minutes and may not be distinguished with the 
capacities of current instruments. The time delay between the NIR and X-ray peaks is also 
within 30 minutes. We therefore expect a good correlation between the NIR and X-ray flux 
densities. The right panel shows the corresponding correlation between ex and the X-ray 
spectral index ax- This correlation gives a correlation between the X-ray flux density and 
spectral index, which can be tested with observations. Interestingly the X-ray flux density 
peaks near a x ~ 0.6 — 0.7, which explains the hardness of X-ray spectra. 

3.4. Implications on the Source Structure 

The above studies clearly show that if the plasma producing the flares does not have 
complicated structure, the observed correlation between the spectral index and flux in the 
NIR band can only be explained by the cooling model with a magnetic field inversely pro- 
portional to the mean energy of the electrons. This is especially true if the total number of 
electrons involved does not change and the electron energy distribution is roughly uniform 
throughout the flare region. Other models either fit the observed correlation marginally or 
predict very weak dependence of the spectral index on the flux in the hard spectral states, 
which is inconsistent with observations. 

On the other hand, the flare region may have complicated source structure. The amount 
of electrons involved in producing the soft state emission can be quite different from that 
producing the hard state emission. If the electron distribution can be approximated as a 
Maxwellian spectrum, as is likely true, the observed correlation can then be used to determine 
the source structure. As shown in section § 3.1, a only depends on B'Jq, and F v also depends 
on BnR 3 . For a given magnetic field, the observed correlation can be used to give a relation 
between the temperature 70 and the total number of electrons at such a temperature, which 
may shed light on the instabilities triggering the flares. One may also assume that the 
sum of the electron and magnetic field pressures is a constant in the flaring region. The 
observed correlation can also lead to a measurement of the amount of plasma at different 
temperatures. 

It is probably more natural to use MHD simulations to produce some flaring regions. 
The formalism developed in this paper can then be used to calculate the radiation spectrum. 
By comparing with the observed correlation, one can then determine what kind of active 
regions are likely responsible for the flares. It is clear that the SSC spectra are quite different 
for all these scenarios and simultaneous spectroscopic observations in the NIR and X-ray 
bands can be used to distinguish these models. Current observations do not justify such a 
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comprehensive investigation, which is also beyond the scope of the paper. 

4. Conclusion and Discussion 

Flares from the direction of the supermassive black hole in the Galactic center are 
produced within a few Schwarzschild radii. Studying the details of the flare properties 
plays an important role in understanding the physics of black hole accretion. To infer the 
underlying physical processes with simultaneous spectroscopic observations over a broad 
frequency range, the plasma process of electron acceleration has to be addressed. Given the 
closeness of the flare variation timescale, the dynamical time at the last stable orbit of the 
black hole, and the synchrotron cooling time of electrons producing NIR emission during the 
flares, a time dependent study of the particle acceleration is necessary. Although there are 
already attempts to carry out this study analytically (Becker et al. 2006; Schlickeiser 1984), 
numerical solutions are required when specific events are considered. 

In this paper we study the time dependent solutions of electron heating and cooling by 
a turbulent magnetic field. When the amplitudes of large scale magnetic fields are much 
smaller than the turbulent component, there are four basic model parameters, namely the 
magnetic field and its coherent length, the gas density and the initial distribution of elec- 
trons. Electrons gain energy from the plasma waves and loss energy via radiation, and the 
energy dependences of the heating and cooling rates are quite different. The evolution of 
the electron distribution is determined by these parameters. We show that for the optically 
thin collisionless plasma in Sagittarius A*, synchrotron cooling dominates and relativistic 
Maxwellian distributions usually give a good description of the electron spectra. A formula 
is given to describe the dependence of the temperature evolution on the model parameters, 
which can be coupled with MHD simulations to study the radiation properties of the plasma. 

We demonstrate that it is not trivial to produce the observed correlation between the 
spectral index and flux in the NIR band. Doppler effects and adiabatic processes suggested 
previously for the correlation (Gillessen et al. 2006; Yusef-Zadeh et al. 2006b) can not 
reproduce the observed results unless one has a very unusual electron spectrum. If the 
magnetic field does not change dramatically during the flare, only the cooling model of 
electrons heated instantaneously at the flare onset fits the observed correlation marginally. 
All other models predict very weak dependence of the spectral index on the flux density in 
the relatively hard spectral states. 

We argue that the dynamical processes of electrons heating and cooling give the most 
reasonable explanation for the observations. To fit the correlation, the magnetic field needs 
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to be anti-correlated with the electron temperature to reduce the flux density in the hard 
spectral states, bringing the above theoretical predictions more in line with the observations. 
It is not clear what mechanism may be responsible for such an anti-correlation. One pos- 
sibility is that the dependence of the heating rate on the magnetic field is not as strong as 
that of the cooling rate. For example, the heating rate may be determined by the sound 
velocity, which can be much higher than the Alfven velocity. Moreover, the cooling due to 
inverse Comptonization will dominate when the electron temperature and column depth are 
very high, which is certainly the case in the early cooling phase 6 . All these effects can give 
rise to the desired anti-correlation. 

The required anti-correlation between the magnetic field intensity and the electron mean 
energy suggests efficient energy exchanges between the field and electrons. While magnetic 
reconnection can cause transfer of magnetic field energy into electrons (corresponding to the 
heating phase), it is not clear what mechanism are responsible for the increase of magnetic 
field in the cooling phase. If the flare region is in pressure equilibrium with a relatively stable 
background, the magnetic field pressure may increase as the electron pressure decreases due 
to cooling. 

We note that the magnetic field pressure increases from 1.27 to 1.27 x 10 4 ergs cm -3 and 
the electron pressure decreases from 6.1 x 10 5 to 6.1 x 10 3 ergs cm -3 in the cooling phase of 
the flare discussed in § 3.2. For the heating phase, the magnetic field pressure decreases from 
3.2 x 10 3 to 0.32 ergs cm -3 , while the electron pressure increases from 2.1 to 2.1 x 10 2 ergs 
cm -3 . Although the pressures are far from equipartition in the initial and final states, the 
observed correlation is produced when the two are evolving toward an equipartition. These 
suggest that the flares are triggered by a process, which drives the energy densities of the 
magnetic field and electrons far from equipartition. The relaxation of the system toward 
energy equipartition produces the observed correlation. It is therefore critical to produce 
and/or identify features far from energy equipartition in MHD simulations to uncover the 
physical processes driving the flares. 

It is also possible that the flare region has complicated structure. The results of the elec- 
tron spectral evolution presented here, which are also applicable to the relatively quiescent- 
state of the accretion flow, then need to be coupled with MHD simulations to study the 
radiation characteristics of the accretion flow. The observed correlation between spectral 



We note that the X-ray luminosity in the early cooling phase is higher than the synchrotron luminosity, 
suggesting a Compton catastrophe. This is due to the choice of the model parameters to embrace a large 
dynamical range. Including the cooling effect due to inverse Comptonization will make the electron temper- 
ature decrease rapidly to a value, where the Comptonization luminosity becomes lower than the synchrotron 
luminosity. 
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index and flux suggests that, if the magnetic field does not change dramatically during the 
flare and through the flare region, more electrons are involved in producing the softer spec- 
tral state flux. Regions with small hot cores flashing in a extended lower temperature active 
envelope in MHD simulations may be identified with the flares. 

This work was carried out under the auspices of the National Nuclear Security Ad- 
ministration of the U.S. Department of Energy at Los Alamos National Laboratory under 
Contract No. DE-AC52-06NA25396. This research was partially supported by NSF grant 
ATM-0312344, NASA grants NAG5-12111, NAG5 11918-1 (at Stanford), and NSF grant 
PHY99-07949 (at KITP at UCSB). SL would like to thank Justin R. Pelzer for developing 
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Log 10 [y] Log 10 [y] 

Fig. 1. — Left: Evolution of the electron distribution function Ni^y). The initial and final 
temperatures (in unit of m e c 2 ) are 10 and 200, respectively. The time steps are 4(i/14) 2 r ac , 
where % is an integer ranging from to 14. The arrow indicates the direction of the evolution. 
The thick line is for i = 7 corresponding to the time step r ac . We note that the electron 
distribution is already very similar to the steady-state spectrum after r ac . Right: Same as the 
left panel except the initial and final temperatures exchanged. Note the pileup of electrons 
during the first few time steps due to the SSC cooling. 
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Fig. 2. — Left: Correlation between N and (3 = dlniV/dm7 at 7 = 50. The solid line is for 
the heating phase. The dashed line is for the cooling phase. The rising of the spectral index 
to above 2 is caused by the pileup of electrons in the early cooling phase. The correlation is 
quite different for the two phases. Right: Same as the left panel but for 7 = 100. 
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Fig. 3. — Left: Evolution of the mean energy of the electrons in the heating phase. The 
dashed line indicates the theoretical result given by equation (6) with a = 1.6. Right: Same 
as the left panel but for the cooling phase with a = 1.0. 
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Fig. 4. — Left: The ratio of the mean energy of the electrons to the theoretical value as 
a function of time for several initial and final temperatures in the heating phase. The 
solid lines have an initial temperature of 10 and final temperatures of 200, 500, 1000. The 
dashed lines have a final temperature of 200 and initial temperatures of 20 50, and 100. The 
relative error (the difference between the numerical and theoretical results divided by the 
theoretical values) increases with the increase of the ratio of the final to initial temperature 
(the dynamical range). Right: Same as the left panel but for the cooling phase. The line 
with the greatest relative error corresponds to a cooling from 1000 to 10. The other lines 
have an initial temperature of 200 and final temperatures of 10, 20, 50, and 100. The relative 
error also increases with the dynamical range. 
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Fig. 5. — Evolution of the normalized synchrotron flux density spectrum e v for the two runs 
shown in Figure 1. The time steps are the same as in Figure 1 and the magnetic field B = 100 
G. The left and right panels correspond to the heating and cooling phases, respectively. The 
evolution of e u is less dramatic than that of N(j) due to the dominance of synchrotron 
emission by more energetic electrons. The spectra are very similar to thermal synchrotron 
spectra. 
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Fig. 6. — Correlation between the normalized flux density e v and the spectral index a at 
v — 1.4 x 10 14 Hz (left) and 1.4 x 10 13 Hz (right) for the two runs in Figure 5. The solid and 
dashed lines are for the heating and cooling phases, respectively. The dotted lines give the 
observational results. The cooling phase correlation fits the observations marginally, while 
the heating phase correlation is too flat compared with observations. See text for details. 



- 22 - 



t 1 1 1 1 1 1 1 1 1 1 1 [ 1 1 1 1 1 1 1 1 1 — | 1 1 1 1 1 1 1 1 r 




— 1_\ 1 1 1 1 1 1 1 1 1 1 1 I ' L- ' I 1 1 1 1 1 1 1 — I 1 1 1 1 1 1 1 1 1 — 

0.2 0.4 0.6 1 2 3 4 



Fig. 7. — Left: Correlation between e u and a for electrons with a relativistic Maxwellian 
distribution. The dotted lines indicate the observed results, which clearly lie below the 
theoretical curve. Right: Correlations between xm^u (solid line), x 2 M e u /50 (dashed line) and 
a, which are more consistent with observations. See text for details. 
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Fig. 8. — Left: Correlation between the flux e u B and spectral index a at 1.4 x 10 14 Hz for the 
heating (solid line) and cooling (dashed line) phases, where the magnetic field B is chosen 
to be inversely proportional to the mean energy of the electrons < 7 >. The initial magnetic 
field, the initial and final temperatures are 4 Gauss, 2000 and 20 and 200 Gauss, 10 and 1000 
for the cooling and heating phases, respectively. The dotted lines indicate the observational 
results. The scale of the heating phase flux is indicated on the upper axis. Right: Time 
evolution of the magnetic field for the heating (solid line; upper scale) and cooling (dashed 
line; lower scale) phases. The regions most relevant to flare (0 to 3 minutes for the heating 
phase and 133 to 163 minutes for the cooling phase) observations are shaded. 
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Fig. 9. — Same as the right panel of Figure 8 but for the flux (left) and spectral index (right) 
at 1.4 x 10 14 Hz. The scale of the heating phase flux is indicated on the right axis. Solid 
circles correspond to the NIR spectra in Figure 10. See text for details. 
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Fig. 10. — Evolution of the SSC spectrum during the cooling (left) and heating (right) 
phases shown in Figure 8. Besides the initial and steady-state spectra, we show the spectral 
evolution for the cooling and heating phases when the correlation is produced. These spectra 
are indicated by the solid circles in Figure 9 and Figure 11. 
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Fig. 11. — Left: Light curve of the flux and spectral index at 0.5 x 10 18 Hz (thin lines, 
corresponding to ~ 2.1 keV) and 2.0 x 10 18 Hz (thick lines) during the cooling (dashed 
lines) and heating (solid lines) phases. The scales for the heating phase are indicated on 
the right and top axes. Right: The corresponding correlation between the X-ray flux and 
spectral index. The scale of the heating phase is indicated on the top axis. The solid circles 
correspond to the SSC spectra in Figure 10. 



